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Abstract — This paper considers linear precoding for constant 
channel-coefficient i^-User MIMO Gaussian Interference Chan- 
nel (MIMO GIC) where each transmitter-i (Tx-i), requires to 
send di independent complex symbols per channel use that take 
values from fixed finite constellations with uniform distribution, 
to receiver-i (Rx-i) for i — 1, 2, • • • , K. We define the maximum 
rate achieved by Tx-i using any linear precoder, when the 
interference channel-coefficients are zero, as the signal to noise 
ratio (SNR) tends to infinity to be the Constellation Constrained 
Saturation Capacity (CCSC) for Tx-i. In this paper, we derive 
a high SNR approximation for the rate achieved by Tx-i when 
interference is treated as noise, which is given by the mutual 
information between Tx-i and Rx-i, denoted by I[Xi;Yi\ where, 
Xi denotes the symbols generated at Tx-i before precoding and 
Yi denotes the symbols received at the antennas of Rx-i. Based 
on this high SNR approximation, we derive a set of necessary 
and sufficient conditions on the precoders under which I[Xi\ Yi] 
tends to CCSC for Tx-i. Interestingly, the precoders designed 
for interference alignment satisfies these necessary and sufficient 
conditions. However, finding precoders that align interference is 
known to be NP-hard in general whereas, the precoders that 
satisfy the derived necessary and sufficient conditions are easy 
to find for any given channel-coefficients. Precoding for the case 
of single antenna at all the nodes correspond to rotation of 
constellation at the transmitters. As a corollary for the single 
antenna case we have the result that, using appropriate rotation 
of constellation at the transmitters, as SNR tends to infinity, 
7[X{;ii] tends to CCSC for Tx-i which is independent of the 
channel-coefficients. This result is in contrast to the Gaussian 
alphabet case where I[Xi\Yi] tends to a value determined by 
the channel-coefficients, as SNR tends to infinity. Further, we 
propose a gradient-ascent based algorithm to optimize the sum- 
rate achieved by precoding with finite constellation inputs and 
treating interference as noise. Simulation study for a 3-user 
MIMO GIC with two antennas at each node with dk — 1 
and QPSK inputs, for all i, shows an improvement of 1.07 
bits/sec/Hz in the ergodic sum-rate using the precoders obtained 
from the proposed algorithm over the precoders that achieve IA, 
at SNR = -2 dB. 

I. Introduction 

Interference Alignment (IA) has been a focus of intense 
research on Gaussian interference channels (GICs) in the 
recent years on the account of the capacity of interference 
channels being unknown in general and the potential of IA to 
get close to the sum-capacity for a broad class of interference 
channels as the signal to noise ratio (SNR) tends to infinity. 
The scaling of the sum-capacity with log SNR is known as 
the sum degrees of freedom (DoF) of the GIC [1|. The sum- 
DoF is known to be K/2, with probability 1, for the K-user 



GIC [1 1 when all the nodes are equipped with a single antenna 
and time-varying channel gains are assumed. The result was 
proved by using linear precoding at the transmitters over an 
arbitrarily large number of symbol extensions and zero forcing 
at the receivers. In a later work J2), the sum-DoF is shown to 
be K/2 with probability 1 even in the case of constant channel- 
coefficients (that are drawn from a continuous distribution), 
i.e., the channel gains do not vary with time, using a non- 
linear IA technique. When all the nodes are equipped with M 
antennas, the sum-DoF was shown to be 3M/2 for the 3-user 
GIC [1|. This result was proved by linear precoding over the 
transmit antennas without the use of symbol extensions and 
holds true even in the case of constant channel coefficients. 
Later [3 1, with the assumption of constant channel coefficients 
and using a non-linear IA technique, the sum-DoF when 
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gcd(M, N) denotes the greatest common divisor of M and N, 
with M being the number of antennas at each transmitter and 
N being the number of receive antennas at each receiver. All 
the works cited above assumed full channel state information 
at all the transmitters (CSIT) and receivers (CSIR). The notion 
of sum-DoF involves scaling of sum-rate as log SNR at 
high SNR and therefore, Gaussian input alphabets or lattice 
codes are always used in the study of sum-DoF. However, in 
all practical scenarios finite constellations like M-QAM and 
M-PSK are used at the inputs. With the constraint of finite 
constellation inputs, it is not known whether IA is optimal in 
some sense. 

Linear precoding for optimizing the mutual information 
between the input and the output has been studied for the 
single user MIMO channel with finite constellation inputs in 
[4|- |6|. Constellation rotation for optimizing the sum-capacity 
for SISO Multiple Access Channel (MAC) with finite constel- 
lation inputs has been examined in |7| and linear precoding for 
weighted sum-rate maximization in MIMO MAC with finite 
constellation inputs has been studied in (8). Note that linear 
precoding for the SISO MAC corresponds to constellation 
rotation at the transmitter. 

Recently, there has been some progress on the analysis of 
finite constellation effects in 2-user SISO GIC l9l. iflOl. In ll9l. 
constellation rotation was found to increase the constellation 
constrained sum-capacity of 2-user SISO Gaussian strong 
interference channel ifTTl . Ifl2l . and in IflOl . a metric to find 
the optimum angle of rotation was proposed. In this paper, 
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we examine achievable rate-tuples with linear precoding for 
if -user MIMO Gaussian Interference Channel (GIC) with 
finite constellation inputs. Specifically, we treat interference 
as noise, i.e., each transmitter reveals its codebook only to its 
intended receiver. The maximum rate achievable under such 
a circumstance for transmitter-i (Tx-i) is given by mutual 
information between the input generated by Tx-i and the 
output at receiver Rx-i. The channel conditions and values of 
SNR under which the decoding scheme of treating interference 
as noise with Gaussian alphabet inputs is sum-capacity optimal 
was found for the 2-user SISO GIC in (B)- Q3), for the K- 
user SISO GIC in QU and for the 2-user MIMO GIC in 07), 
|18|. For given values of channel gains, with Gaussian input 
alphabets, as the SNR tends to infinity, treating interference 
as noise is not sum-capacity optimal |19|. 

With the constraint of finite constellation inputs, it is not 
clear whether treating interference as noise is optimal in 
some sense. First, we need to define a notion of optimality 
under the constraint of fixed finite constellation inputs and 
then analyse decoding and transmit schemes with that notion 
of optimality. Consider a scenario where each transmitter-i 
(Tx-i), requires to send d, independent complex symbols per 
channel use that take values from fixed finite constellations 
with uniform distribution, for i — 1, 2, ••• , if, to receiver- i 
(Rx-i). Throughout this paper, we assume that none of the 
direct channel gains are zero. For a if -user MIMO GIC with 
finite constellation inputs, as a measure of optimality of linear 
precoding in the high SNR sense, we introduce the notion of 
Constellation Constrained Saturation Capacity (CCSC) which 
is defined as follows. 

Definition 1: The maximum rate achieved by Tx-i as SNR 
tends to infinity, using any linear precoder, when the interfer- 
ence channel-coefficients are zero is termed as the Constella- 
tion Constrained Saturation Capacity (CCSC) for Tx-i. 
For the ease of exposition, throughout the paper, we assume 
that the constellations used for the symbols are all the same at 
all the transmitters, and is of cardinality M. Hence, the CCSC 
for Tx-i is given by logiM di . 

In this paper, with the assumption of constant if-user 
MIMO GIC with full global knowledge of channels gains, 
and finite constellation inputs, we derive a set of necessary 
and sufficient conditions on the precoders under which treating 
interference as noise at Rx-i will achieve a rate for Tx-i that 
tends to CCSC for Tx-i, for all i, as SNR tends to infinity. 
Precoders satisfying these necessary and sufficient conditions 
exist for all direct and cross channel gains, and are termed 
as CCSC optimal precoders. Hence, in the case of finite 
constellation inputs with the use of appropriate precoders, the 
rate tuples obtained by treating interference as noise tend to 
values that are independent of the channel gains. For a K- 
user SISO GIC, this result is in contrast with the Gaussian 
input alphabet case where the rate tuples obtained by treating 
interference as noise tend to values dictated by the channel 
gains, as the SNR tends to infinity. Interestingly, the precoders 
that achieve IA, if feasible, are also CCSC optimal precoders. 
However, finding precoders that align interference is known 



to be NP-hard ll20l in general whereas, the precoders that 
satisfy the derived necessary and sufficient conditions are easy 
to find for any given channel-coefficients. Since finite SNR is 
of more practical interest, we propose a gradient-ascent based 
algorithm to optimize the precoders for the sum-rate achieved 
by treating interference as noise. The contributions of the paper 
are summarized below. 

• For a constant if -user MIMO GIC using finite 
constellation inputs with precoding, a high SNR 
approximation for the rate tuples achieved by treating 
interference as noise at the receivers is derived (see 
Theorem Q] in Section [TTTb. Based on this approximation, 
we derive a set of necessary and sufficient conditions 
under which the precoders are CCSC optimal (see 
Theorem [2] Section UTTb . These conditions are satisfied 
with probability 1 when the entries of the precoders are 
chosen from any continuous distribution. It is observed 
that the precoders that achieve IA, if feasible, are CCSC 
optimal. 

• For the finite SNR case, we propose a gradient-ascent 
based algorithm to improve the sum-rate achieved by 
treating interference as noise using finite constellation 
inputs with precoding. Simulation studies indicate con- 
siderable improvement in the ergodic sum-rate with pre- 
coders obtained using the proposed algorithm for a 3-user 
MIMO GIC with 2 antennas at all the nodes, d, = 1, for 
all i, and QPSK inputs, at low and moderate SNRs over 
that obtained using the IA solution of (TJ. 

The paper is organized as follows. The system model is 
formally introduced in Section [TTJ In Section [Till a set of 
necessary and sufficient conditions for CCSC optimal pre- 
coders is derived. In Section IIV1 our gradient-ascent based 
algorithm for optimizing the sum-rate using precoders and 
treating interference as noise at the receivers is given, and 
a simulation result comparing its performance with respect 
to that of precoders satisfying IA is presented. Section [V] 
concludes the paper. 

Notations: For a random variable X which takes value from 
the set X, we assume some ordering of its elements and use x l 
to represent the i-th element of X. Realization of the random 
variable X is denoted as x. The notation diag{V\, Va, • • ■ , V n ) 
denotes a block diagonal matrix formed by the matrices V{, 
i = 1,2, • • • , n. The i th coordinate of a complex vector X 
is denoted by X(i). The 2-norm of a complex vector X is 
denoted by The cardinality of a set X is denoted by \X\. 
For a complex number a, 3?{a} and 3{a} denote the real and 
imaginary parts of a respectively. For two complex numbers 
a and b, the notation a > b denotes that |5i{a}| > |5ft{£>}| and 
> The notation represents the zero vector 

whose size will be clear from the context. All the logarithms 
in the paper are to the base 2. 

II. System Model 

Tx-i intends to communicate with Rx-i, for i = 1, 2, • ■ • , if , 
as shown in Fig.Q] Without loss of generality, let P denote the 
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Fig. 1. System Model. 



power constraint at all the transmitters. The signal received at 
Rx-j is given by 



Yj =J2VPH ij V i Xi 



•JVi 



where, Hij denotes the constant channel matrix from Tx-i 
to Rx-j, Vi denotes the precoder at Tx-i, Xi denotes the 
complex symbol vector generated at Tx-i, Nj denotes the 
noise random vector whose coordinates represent independent 
and identically distributed zero mean unit variance circularly 
symmetric complex Gaussian random variables. The sizes of 
the matrices Hy, Vi, Xi, and Nj are given by n r . x n ti , 
n ti x di, di x 1, and n r . x 1 respectively, where n r , and 
n ti denote the number of receive and transmit antennas at 
Rx-j and Tx-i respectively, and di denotes the number of 
independent complex symbols per channel use that Tx-i wants 
to transmit to Rx-i. These complex symbols are assumed to 
take values from finite constellations with uniform distribution 
over its elements. For simplicity of exposition, we assume 
that the finite constellations used are all the same at all the 
transmitters and are of cardinality M . The results of this paper 
apply with simple modifications when this is not the case. The 
finite constellation used is denoted by S and is of unit power. 

III. CCSC Optimal Precoders 

In this section, we shall derive a set of necessary and 
sufficient on the precoders for CCSC optimality which is 
taken to be a measure of optimality for linear precoding in 
the high SNR regime for the finite constellation input case. 
Rate achievable for Tx-i by treating interference as noise at 
Rx-i, for all i, is given by Ri < I[Xi\Yi\. Our focus will 
be on the boundary point given by I[Xi]Yi\, for all i. Let 
X = [Xi X2 ■ ■ ■ Xk] t ■ The effective channel matrix from all 
the transmitters to Rx-i is given by Hi = [Hu H2i ■ ■ ■ Hki\- 
Define 

V = diagiVuVz,- ■■ ,V K ). 
Using the chain rule for mutual information lETl . 



With uniform distribution assumed over the elements of the 
constellation, the expression for I\X\,X2,--- ,XK\Yi\ is 
given by (O (at the top of the next page) which is derived 
in a similar way as in iflOl . Define the matrix 



A 



HiV(x 



fei 



The following theorem gives a high SNR approximation for 
I[Xi, X 2 , • • • , Xk] Yj\. 

Theorem 1: At high P, the mutual information 
I[Xi, X2, • • • , Xk; Yi] can be approximated by 
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Proof: Define the sets 

A kl = {k 2 ^ h I A klM = 0} 
A klM = {I I A klM {l) = 0} . 



(10) 
(11) 



The expression in is now re-written as in (0) and ©. 
At high values of P, the following inequality holds good for 
k 2 hi, k 2 £ A kl , and for alU i A klM . 



\$t{VPA klM (l)}\ » 



A and \Q{^PA k ^(l)}\ » -|. 



(12) 

Now, note that for k 2 ^ fa, k 2 £ A kl , and for all I £ A klM , 
^/PA klM (l) » Hj(l) when \n ]R (l)\ < 3 and \n 3l (l)\ < 3 
where, rijji(l) and n,ji{l) represent the real and imaginary 
parts of rij(l). The value of rij(l) becomes comparable to 
y/PA klM (l) only if either \n jR (l)\ > 3 or |n i7 (/)| > 3. How- 
ever, the probability of such an event occurring is extremely 
small because the variances of the Gaussian random variables 
rijji(l) and are equal to |. Hence, the contribution 

of such an event to the integral in (O is very small that it 
can be neglected at high P. Therefore, on account of (fl2l) . 
the approximations in © and (0 are valid. The equation (HJ 
follows from the fact that probability distribution integrates to 
1, (O follows from ( fTTT i. and the proposed approximation in 
(fl} is obtained directly by re-writing ([9). ■ 
Remark 1: Note that we cannot straightforwardly argue that 
at high powers P, the noise N in (O can be neglected. This 
is because the value of ||A^|| can be of the order of ||^4' Cl,fe2 1|. 
The proof uses the fact that the probability of such an event is 
very small and hence, can be neglected. The greater the power, 
the better is the approximation. A similar approximation was 
developed for the SISO case in IflOl . The approximation given 
in Theorem Q] is a generalization of the approximation derived 
in d. 
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matrix from all the transmitters, with the exclusion of Tx-i, to 
Rx-i is given by 



I[Xi;Yi\= I[X!,X 2 ,--- ,X K ;Yi\ - J[Xi,AV-- ,X K ;Yi\X t }. 



Hf — [Hn H12 ■ ■ ■ Hi-i^i H i+ i :i ■ ■ ■ PiKi 



I(X 1 ,X 2 ,--- ,X K -Yi) = logM Y -i=^ di 

M E i = l d i 



E E » 



log ^ exp - 

k 2 =0 



■ logM^*=l di 



logM^i=l di 
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(5) 
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it j — r 
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Define the matrices 



B i *' i > = HfVf(x i }-xj), 

for ii, i 2 = 0, 1, • • • , M^i** dj - 1. Similar to (Q), we have 
the following approximation for I[Xi,X2, ■ ■ ■ ,Xx;Yi\Xi] at 
high P. 



log M J =1 

A/ E J^ d J-l 

21 —0 
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Now, define the set 

S J1 = {i 2 7^ ii | B h ' i2 = 0} . 



(15) 



The following theorem gives a set of necessary and sufficient 
conditions under which the above approximation tends to 
Hence, a high SNR approximation for /[X^l^] is given by log M di as P tends to infinity and hence, gives a set of 



necessary and sufficient conditions under which the precoders 
are CCSC optimal. 

Theorem 2: The approximation for I[Xi\Yi] given in ( TBI 
tends to log M di as P tends to infinity iff 

K 

HuViixf 1 - xf 2 ) + J2 HkjVkixl* 1 - x™ 2 ) £ 0, (16) 

k = l 

V Pn + Pi2, V p k l, Vp fc2 

where, p kl ,p k 2 = 0, 1, • • ■ ,M d < -1. 

Proof: The summation-term of the second term in ( TBI ) is 
re-written as 



M^l".-l 
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l + |.A fc i| + E e 
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Similarly, the summation-term of the last term in ( TT4i > tends 

to 53ff=o ^ (l + I) as P tends to infinity. Hence, as 
P tends to infinity, (fl"4l > tends to 



m2-;=i d - 



J2 log(l + \A^\) 



(17) 
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52 ^(l + IB' 1 !). 



Now, define the sets 



C P11 '"' * PK1 = { (P12.P22, ■ ■ ■ ,PK2) ¥= (P11,P21, ■ ■ ■ ,PXl) I 
K 

HiiViitf* - <* 2 ) + J2 H kjV h {x v k " - < fc2 ) ^ }, 

k^i 

J5P11,"' >Pi-l,l>Pi + l,lr" .Pifl 

= { (P12, • • • ,Pi-l,2,Pi+l,2, ■ ■ ■ ,PK2) 

+ (Pll, ■ ■ ■ ,P»-l,l,Pi+l,l> • ' ' iPXl) I 

J2 H kjVk(x p k "-xl*i)^0} 

kj^i 
k=l 

where, Pn,pi2 = 0, 1, • • ■ , M dl — 1. Observe that the set of all 
qpii, ■■■ ,pki }j as a one-one correspondence with the set of all 
A k \ and the set of all ppu>- .Pi-i.i.w+i.i.- ,pki has a one- 
one correspondence with the set of all B n . Hence, ( fTTT i can 
be re-written as 



log M" 



1 



M a l-1 M a Jf-l 



£ ... J2 log(l + \CV 



ll.'" ,PJC 1 
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P11=0 P.Ki=0 
M d l-1 A/ di -l -1 A/ di + l -1 M d K -\ 

. V V V ... V (is) 

P11=0 Pi-l,l=0 Pi+l,l=0 PfCl=° 

iog(l + [DPii.-.«-i.i.p»+i,i.-- «i|) . 



Now, the set C Pl1 ' "' ' P/fl can be written as a disjoint union of 
two sets, i.e., 

QPll,--- ,PK1 _ QPll,--- ,PK1 y (JPH'"' 'P-^ 1 (JC)) 

where, Cf 11 '"' pjn and Cf 11 '" ,p,cl are defined in © and 



respectively (given at the top of the next page). 
The set C^ U '"' ,J ' K1 can be re-defined as in (f2Tb . Since 



Pkl „Pk2\ 



is independent of the indices 



PiX and pj2, the set (jP 11 ' "^^ 1 j s the same for all 
p a = 0, l,--- ,M di - 1. Now, note that the set of all 

£f>n, ,p»-i,i, ,p.+i,i, ,pki nas a one . one correspondence 

with the set of all £>ph>- ,Pi-i,i,Pi+i,i,- ,pjti w hi c h follows 
from the definitions of the respective sets. Hence, ( fT~8b can be 
re-written as 



log M di 

_ uELi. E " E E E ■■■ E 

pil=0 pi_i i i=0 Pi, 1=0 Pi+i,i=0 PK1=0 

Zoc, (1 + |X> Pl1 '- .P»-i,i>Pi+i,i.- ,PK1 | + [Cf 11 '"' 1) 

^ A/ d l-l A/ d '-l -1 Af d i+1 -1 A/ d /f-l 

M P11=0 Pi _i,!=0 p I + i,i=0 P/fi=0 

Zog(l + |X)Pii.--.Pi-i,i.P.+i,i.-".PKi|) 



A/ d l-l AO-l-UA-lJiri+l-l M d Jf_l 

E E E E ••• E 

pn=0 P»-l,l=0 Pi, 1=0 Pi+l,l=0 PK1=0 



log(l + |DPn.-.Pi-i,i.Pi+i,i.-.PKi| + | C Piir-,PKi|) 

M d l-1 A/ d »-l -1 A/ d i -1 A/ d >+l -1 M d K-l 

E •■■ E E E ••• E 

P11=0 P»-1,1=0 Pi, 1=0 p i+ i,i=0 Pki=0 

log {I + [2>Pii,-,Pi-i,i,Pi+i,i,-,Piri[), 



Clearly, if \C, 



Pll,— ,PK1 I 



> 1 for some (pn, ■ ■ ■ ,Pki) then, the 



second term in the above equation is strictly greater than the 
last term and hence, as P tends to infinity, I[Xi \ Yj\ tends to a 
value that is strictly less than log M d \ If jcf 11 '"' ,PK1 1 = for 
all pn then, the second term in the above equation is equal 
to the last term and hence, /[JQ;!^] tends to a value equal 
to log M d \ Thus, I[Xi;Yi] tends to log M dl as P tends to 
infinity iff ( fT6] l is satisfied. ■ 
Remark 2: The result of Theorem |2] means that the rate 
achieved by treating interference as noise at high P tends to 
CCSC for Tx-i iff, in the absence of the Gaussian noise, two 
different symbol vectors x\ !l and x^ 2 sent by Tx-i should not 



c pii.- -.pki = { ( pi2:P22i ... tPK2 ) £ (pu.pa,... ,p K1 ) forp i2 = Pil | HuViix^ 1 - x***) + J2 H k] V k (x P k kl -x P k k2 )=0}, 

k = l 

K 

=>. CfU'- ^A" 1 = { (pi 2 ,p 22 ,-- • , P K2) tL (P11,P21," - ,PKl) forp i2 = Pil | ^^^K" ~ ^ ) = fi } 



fc=l 



C *11.-.*K1 = { (pi2:P22i ... >pK2 ) ^ ( Pll ,p 21 ,... ,p K1 ) forp i2 # Pll | ffiiVUxf" -zf^ + ^^VU^" ~x P k " 2 ) = } 



(201 



(21) 



(22) 



map to the same symbol vector at Rx-i for any data symbol 
transmitted by the interfering transmitters. 

Remark 3: For a given value of channel gains with none 
of the direct channel gains being 0, when the entries of the 
precoders are chosen from any continuous distribution (say, 
standard normal distribution) the probability of the event 



HaVi(xf l - xf 2 ) + J2 H kjV k (x P k kl 
fe=l 



o, 



to occur for any pn ^ p i2 and for any (pki, Pk2) is zero. 
By appropriate scaling of the precoders thus obtained, with 
probability 1, we have CCSC optimal precoders. 

Remark 4: Interference alignment, if feasible ll22l for the 
given values of nt t , n n , and di involves finding precoders 
such that the signal sub-space at Rx-i, generated by [HuVi], 
is linearly independent of the interference sub-space, generated 
by [H u Vx ■■■ II, i ,i; : //,. ; ••• H K .iV K }, and the 
matrix [HuVi] is full-rank Q']. The CCSC optimality condition 
in ( fT6b can be rewritten as 




— Tx-1 with QPSK 
■a-Tx-2 with QPSK 
-•-Tx-3 with QPSK 
--Tx-1 with Gaussian Alphabet 
+-Tx-2 with Gaussian Alphabet 
-•-Tx-3 with Gaussian Alphabet 
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Hi-i^Vi-i H i+ i }i V i+ i ■■■ H Kii V K ] 



PKI 
K 



,:pk2 



^ 0, V pil ^ pi- 



Since, with precoders that achieve IA, the signal sub-space 
at Rx-i is linearly independent of the interference sub-space 
and [HuVi] is full-rank, the above condition is satisfied for 
all i. Hence, IA precoders are also CCSC optimal precoders. 
However, in general, finding such precoders are NP-hard [20 1 
whereas finding CCSC optimal precoders are easy to find as 
explained in the previous remark. 

Example 1: Consider a MIMO GIC with K = 3, n ti = 
n ri = 1, dj = 1 for all i, and the finite constellation used is 
QPSK. The channel matrix and the precoders are given by 

"-0.9 + 0.4i -l.7-l.40i 1.5 + 5.0i 
H= 2.6-0.9i -0.9-2.8i 0.04 + 0.88t 
-2.9- 5. 2i -10.2 + 0.7i -0.5 + 2.4i] 

V, = l,V 2 = ,V 3 = 1 

where, the matrix element represents the channel gain 

from Rx-i to Tx-j. The mutual information I[Xi \ Yi\ evaluated 



Fig. 2. Rates in bits/sec/Hz vs P in dB for Example [T] 



using Monte-Carlo simulation is plotted for QPSK inputs and 
Gaussian inputs, for all i, in Fig. [2] The chosen precoders 
satisfy ( TToT ) and hence, I\Xi\ Yi\ saturates to 2 bits/sec/Hz for 
all i, as P tends to infinity in the QPSK case whereas, for 
the Gaussian alphabet case, the saturation rate is determined 
by the channel gains. The saturation value of I[Xi] Yi] in 
the Gaussian alphabet case is given by log ^1 + ^ ^'|L fc | 2 ) 
which evaluates to 0.04, 1.02, and 0.06 bits/sec/Hz for Tx-1, 
Tx-2 and Tx-3 respectively in this example. 

Example 2: Consider a MIMO GIC with K — 3, n u = 
n ri = 2, di = 1 for all i, and the finite constellation used is 
QPSK. Let the effective matrix from all the transmitters to all 
the receivers be given by 





'H n 


H21 


H31 


H = 




H22 


H32 




Hv3 


H23 


H33 



where, Hij is the 2x2 channel matrix from Tx-i 
to Rx-j. The precoders and the channel matrices are 
given in d23l ) and d24b respectively (given at the top 
of the next page). The QPSK points are given by 

Now, (O 



75 + y/2 1 ' 



1 
V2 



V2 L 



1 

V2 



V2 l 



1 

V2 



72' 



H = 



0.66 + 0.74i 
0.13 + 0.99i 

0.5756 - 0.0565i 
0.1610 + 0.3766i 
-1.1533 - 0.1280i 
-1.7763 - 0.3748i 
-1.7082 - 0.4948i 
-1.3014 - 0.5614i 



.9883 + 0.1524i 
4538 + 0.891H 

0.7524 - 0.1375i 
-0.0010 + 0.2005i 
-0.6361 + 1.4658i 

0.5341 + 0.0966i 
-0.6101 - 0.4739i 

1.2515 + 0.3414i 



0.1697 - 0.1069i 
0.8758 - 0.0689i 
-1.3069 + 0.1090i 
-0.9491 + 0.8074i 
-0.2226 - 4.2486i 
0.4242 + 0.0202i 



7044 + 0.7098i 
1603 + 0.987H J ' 

0.0124 - 0.2002i 
-0.1285 + 0.0605i 
0.0427 + 0.2488i 
-1.0773 - 1.7202i 
-0.8216 + 0.4808i 
0.0138 - 0.8740i 






-0.0028 + 0.2215i 
0.9616 - 1.2130i 
0.9572 + 1.8870i 
0.3393 - 1.345H 






-1.0597 - 0.2708i 
-0.6077 + 0.6970i 
-1.4428 - 1.4353i 
0.9498 - 1.0932i 



(23) 



(24) 




1( R Monte Carlo 


R Approximation 


-$-R 2 Monte Carlo 


-a- r 2 Approximation 


_I_R 3 Monte Carlo 


t R^ Approximation 



max/(Vi, • • • ,Vk) = max^lfX,; Yi\ with Tr^V") < 1. 
j=i 

This is a non-concave problem in general and difficult to 
solve. Hence, we propose a gradient-ascent based algorithm to 
improve the sum-rate starting from some random initialization 
of precoders. Define the MMSE matrix at Rx-j by 



E 4 



i[(X--E[X\Y j ]){X-K[X\Y j ]) H ] 
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where, E represents the expectation operator. Define the 
MMSE matrix at Rx-j with the exclusion of Tx-j's signal 
by 

E f = E[(X f - E[X f \Yj - HjjXj\)(Xj - E[X f \Yj - HjjXj]) 11 ] 

The gradient of the sum-rate with respect to the precoder Vi 
given by 



Fig. 3. Rates in bits/sec/Hz vs P in dB for Example ff] 



is not satisfied for i = 1 because 



[tfnVi H21V2 H31V3] x 



(x^ 



■ 0, for 



(x?» 

(4 31 



X 



P12 \ 
1 J 

4 a2 ) 





+ \/2i 












For i = 2,3, ( TTol l is satisfied. The plots of I[Xi; Yi] evaluated 
using Monte-Carlo simulation and I[Xi\Yi] evaluated using 
the high SNR approximation in (fl4l i is shown in Fig. [3] Note 
that i?i saturates to a value strictly less than 2 bits/sec/Hz 
whereas R 2 and i? 3 saturates to 2 bits/sec/Hz thus validating 
Theorem [2] 

IV. Gradient Ascent Based Algorithm for 
Finite-SNR 

In the previous section, we studied the rate achieved for each 
transmitter by treating interference as noise at every receiver 
as P tends to infinity. In this section, we focus on the finite 
SNR case. Specifically, the aim is to maximize the sum-rate 
achieved by treating interference as noise at every receiver 
with respect to the precoders, i.e., 



i=K 

WJ(Vi, • • • , V K ) = Vv, T l X f> Y A 

3=1 

i=K 

^Vy^/^Ij,... ,X K ;Yj]-I[X 1 ,X 2 ,-- ,X K \Y S \X S ] 
3=1 

i=K 

(25) 

3=1 
-K 



3=1 

log e 22 H i:j H f E f it-JjUi-jl+i^U d k -d j i(i-j) 



3=1 



where, d25l l follows from the relation between the gradient of 
mutual information and the MMSE matrix obtained in [4|. The 
matrices 



and 



select the column numbers from Y?/Z?i dk + i to Y?k=i dk °f 



Ej and Y] 



i-i 
fc=i a * 



djl{i + l to £J. =1 dk ~ d^i" 1 - 3) of 



Ef respectively, where 



i > j 
i < j. 



Define V — diag(Vi, V-j, • • • , Vrt). The gradient ascent 

based algorithm for optimizing /(Vi,-- - , Vk) with 

respect to the precoders is given in Algorithm 1. 

During every iteration, whose number is denoted by 

n, all the precoders are updated as given in Line 



represents 



10 of Algorithm 1 where, Vy/|y = y(„ 

••• ,Vv K f\y K=v (p-i) 
and Vvj/L, „(n-ij denotes the gradient Vy t / evaluated at 
A (n-i) ^ p OW6r cons traint for any transmitter Tx-i 



diag[VvJ\ Vi=v (* 



-i)>Vv 3 /| Va=v tn-i), 



is violated then, ; is projected onto the feasible set with 
Tr{v} n) v£ n)H ) < 1 (see Line 12 of Algorithm 1) [4J. The 
condition in Line 15 of the algorithm ensures that there is 
sufficient increase in the objective function. The step size 
t of the algorithm is chosen by back-tracking line search 
with parameters a and /3 whose typical values lie between 
(0.01,0.3) and (0.1,0.8) (23]. The proposed algorithm stops 
when either the number of iterations performed is equal to 
max iterations or — /(™~ 2 ) < g (see Line 5 of 

Algorithm 1), for some fixed e. 

Algorithm 1 Gradient Ascent based Algorithm for improving 
sum-rate 

Initialize Vj = V z (0) with Tr(V~;V; H ) < 1, i = 1, 2, • • • , K, and t = 1. 



for n = 1 to max_iterations do 



Compute = /(V^ 



i?'"- 1 ', and E^- 1} , fori 



if n > 1 and - /f«- 2 ) < e then 

exit for 
end if 

Compute Vy/|y = y(„-l) 

do 

V (n) ^y(n-l) +tVv/| v =vC»-D- 



1,2, 



if Tr(y/"' V: ( " )H ) > 1, for all • 



Compute /<") = /(V" 1 t ™ ) , V 2 l 



t = /3t 

while /(») < /(«"!) + at ||Vy/|y = y( n _i) 1 1* 
i = 1. 
end for 



Similar gradient ascent based algorithms have been pro- 
posed in the past for optimizing rates in single user MIMO 
channels g), Q and MIMO MAC [8| with finite constellation 
inputs and precoding. Like in [4| [5| [8|, the algorithm does 
not assume uniform distribution over the elements of the 
finite constellation. In the following example, we present some 
simulation results with the proposed algorithm. 

Example 3: Consider a MIMO GIC with K = 3, n ti = 
n ri = 2, di = 1 for all i, and the finite constellation used is 
QPSK. 

1) Performance for a fixed channel: The parameters used 
in Algorithm 1 are a = 0.05, /3 = 0.4, e = 0.01 and 
max iterations = 15. The chosen channel matrix is 
given in d26l l. The convergence behaviour of Algorithm 
1 is shown in Fig. for P = -2dB, P = OdB, 
and P = 2dB with IA precoders of HI used as 
initialization. The proposed algorithm terminates well 
before the max iterations number for all P because 
the condition in Line 5 of Algorithm 1 is satisfied. The 
sum-rate improvement obtained on application of the 
proposed algorithm with IA precoder initialization is 
observed to be 2.35 bits/sec/Hz, 2.12 bits/sec/Hz, and 




3 4 5 

Iteration Number 



(a) P - 



-2dB 




4 5 6 7 

Iteration Number 
(b) P = OdB 




5 6 7 
Iteration Number 

(c) P = 2dB 



Fig. 4. Sum-rate increase with every iteration for the constant channel case 
in Example. [3] 

1.5 bits/sec/Hz at P = -2dB, P = OdB, and P = 2dB 
respectively. 

2) Performance averaged over channel realizations: The 
ergodic sum-rate (i.e., sum-rate averaged with the entries 
of the channel matrices being taken from CA^(0, 1)) 
obtained by IA precoders of [ 1 1 and precoders obtained 
from Algorithm 1 with the IA precoders as initialization 
is shown in Fig. [5] The parameters used in Algorithm 1 
are a = 0.05, (3 — 0.4, e = 0.01, max iterations = 8. 
The sum-rate increase obtained using the algorithm is 
1.07 bits/sec/Hz at P = —2dB which then decreases 
as P increases, and at P = lOdB the sum-rate in- 



" 1.1408 - 0.8637i 
0.9331 + 0.3749i 
-0.5206 + 2.72244 
-0.5897 + 0.83224 
1.3814 - 0.4690i 
.-0.3043 - 0.2351i 



0.1954 - 1.5172i 
1.1588 + 0.1093i 
0.0330 + 1.0577i 
-0.9552 + 1.1568i 
-0.2731 - 0.63554 
1.1578 + 0.3676i 



-0.7038 + 0.30644 
1.1412 + 0.05614 
-0.1581 - 0.7055i 
-0.9685 - 0.724H 
-0.8927 + 0.00444 
0.9095 - 0.4187i 



0.3323 - 0.4278i 

0.2708 + 0.53204 
-1.4823 + 0.7144i 

1.1944 - 0.0556i 
-0.1312 + 0.31064 

0.4161 + 1.2482i 



1.3425 - 0.58424 
-1.1639 + 0.3706i 

0.4286 + 0.5008i 
-1.3293 - 0.2445i 
-1.9229 - 0.2353i 

1.2900 - 0.0196i 



0.0766 - 0.75554 " 

0.4311 + 1.47264 
-1.0052 + 0.00514 

0.3540 - 0.76444 
-0.7885 - 0.32294 

0.4838 + 1.45784 . 



crease is just 0.1 bits/sec/Hz. The negligible sum-rate 
improvement at high values of P is expected because IA 
precoders are CCSC optimal precoders which saturate 
the sum-rate at 6 bits/sec/Hz. 




3.5 



— Gradient Ascent Solution 
- - - Interference Alignment solution 
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Fig. 5. Ergodic sum-rate in Example [5] 

In (TJ, for sum-DoF analysis zero-forcing matrices are used 
at every receiver. Here, we have treated interference as noise 
and have not used any zero-forcing matrices as using them 
can only reduce the rate I\Xi \ Yj\ because of data-processing 
inequality l2D . 

V. Conclusion 

The paper discussed linear precoding for i\~-user MIMO 
GIC with finite constellation inputs. We showed that, for 
constant MIMO GIC with finite constellation inputs, CCSC 
for every transmitter can be achieved just by using a naive 
scheme of treating the interference as noise at every receiver, at 
high SNR. This result is in contrast with the Gaussian alphabet 
case where, at high SNR, the scheme that treats interference as 
noise saturates to a value determined by the channel gains for 
the SISO case. A set of necessary and sufficient conditions 
for CCSC optimal precoders were derived. It was observed 
that IA precoders fall under the class of CCSC optimal 
precoders. However, CCSC optimal precoders are easy to 
obtain for any given value of channel gains unlike obtaining 
IA precoders. Finally, a gradient ascent based algorithm was 
proposed to improve the sum-rate achieved starting from any 
random initialization of precoders. An example simulation 
with QPSK input alphabets showed a significant improvement 
in the ergodic sum-rate achieved by the precoders obtained 



from the proposed algorithm compared to the ergodic sum- 
rate obtained from IA precoders, treating interference as noise 
at the receivers, at low and moderate SNRs. 
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